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Abstract 

In this paper we investigate cosmological dynamics on the normal branch of a DGP- 
inspired scenario within a phase space approach where induced gravity is modified in the 
spirit of /(i?) -theories. We apply the dynamical system analysis to achieve the stable 
solutions of the scenario in the normal DGP branch. Firstly, we consider a general form 
of the modified induced gravity and we show that there is a standard de Sitter point in 
phase space of the model. Then we prove that this point is stable attractor only for those 
f{R) functions that account for late-time cosmic speed-up. 
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1 Introduction 



There are many astronomical evidences supporting the idea that our universe is currently 
undergoing a speed-up expansion [1]. Several approaches are proposed in order to explain the 
origin of this novel phenomenon. These approaches can be classified in two main categories: 
models based on the notion of dark energy which modify the matter sector of the gravitational 
field equations and those models that modify the geometric part of the field equations generally 
dubbed as dark geometry in literature [2,3]. From a relatively different viewpoint (but in the 
spirit of dark geometry proposal), the braneworld model proposed by Dvali, Gabadadze and 
Porrati (DGP) [4] explains the late-time cosmic speed-up phase in its self-accelerating branch 
without recourse to dark energy [5]. However, existence of ghost instabilities in this branch of 
the solutions makes its unfavorable in some senses [6]. Fortunately, it has been revealed recently 
that the normal, ghost-free DGP branch has the potential to explain late-time cosmic speed-up 
if we incorporate possible modification of the induced gravity in the spirit of /(i?)-theories [7]. 
This extension can be considered as a manifestation of the scalar-tensor gravity on the brane. 
Some features of this extension are studied recently [8,9]. 

Within this streamline, in this paper we study the phase space of the normal DGP cosmology 
where induced gravity is modified in the spirit of /(i?)-theories. We apply the dynamical system 
analysis to achieve the stable solutions of the model. To achieve this goal, we firstly consider 
a general form of the modified induced gravity. We obtain fixed points via an autonomous 
dynamical system where the stability of these points depends explicitly on the form of the 
/(i?) function. There are also de Sitter phases, one of which is a stable phase explaining the 
late-time cosmic speed-up. Secondly, in order to determine the stability of critical points and 
for the sake of clarification, we specify the form of f{R) by adopting some cosmologically viable 
models. The phase spaces of these models are analyzed fully and the stabihty of critical points 
are studied with details. 



2.1 The basic equations 

Modified gravity in the form of /(i?) -theories are derived by generalization of the Einstein- 
Hilbert action so that R (the Ricci scalar) is replaced by a generic function f{R) in the action 



where Lm is the matter Lagrangian and = SttG. Varying this action with respect to the 
metric gives 



2 DGP-inspired f{R) gravity 




(1) 




where T^^^ — diag{p, —p, —p, —p) is the stress-energy tensor for standard matter, which is 
assumed to be a perfect fluid and by deflnition /' = ^ . Also Tj£> is the stress-energy tensor 
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of the curvature fluid that is defined as follows 



t!£> = ^g,u [f{R) - Rf] + f -^"^{9a,gfi. - Qapg,.) ■ 



(3) 



By substituting a flat FRW metric into the fleld equations, one achieves the analogue of the 
Priedmann equations as follows [10] 



3fH' = K'pm + [Uf{R) - Rf) - ^Hf 



-2fH = AcVm + R^f" + {R- HR)f" , 



(4) 
(5) 

where a dot marks the differentiation with respect to the cosmic time. In the next step, following 
[9] we suppose that the induced gravity on the DGP brane is modified in the spirit of f{R) 
gravity. The action of this DGP-inspired f{R) gravity is given by 



(6) 



where qab is the five dimensional bulk metric with Ricci scalar TZ, while qab is induced metric 
on the brane with induced Ricci scalar R. The Priedmann equation in the normal branch of 
this scenario is written as [9] 



if) 



3H 



(7) 



where rc — = ^ is the DGP crossover scale with dimension of [length] and marks the IR 
(infra-red) behavior of the DGP model. The Raychaudhuri's equation is written as follows 



^ 2HrJ' 



R^f" +{R- HR)f 
2f 



(8) 



To achieve this equation we have used the continuity equation for p^^^ as 



Rf" 



K'pmRf" 



where the energy density and pressure of the curvature fluid are defined as follows 

1 n 



(9) 



p 



if) 



f{R)-Rf -3Hf'), 



p(^) = \ (2Hf' + f'-l \f{R) - Rf 



(10) 



(11) 



After presentation of the required field equations, we analyze the phase space of the model 
fully to explore cosmological dynamics of this setup. 
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2.2 A dynamical system viewpoint 

The dynamical system approach is a convenient tool to describe dynamics of cosmological 
models in phase space. In this way, we rewrite equation (7) in a dimensionless form as 

. ^rn L (12) 

3H^f' HrJ' GH^f 6H^ Hf ^ ^ 

In the present study, we firstly consider a generic form of the f{R) function, so that one can 
define the dynamical variables independent on the specific form of the f{R) function as follows 
(see for instance Ref. [10]) 

Pm 1_ -_JL nq^ 

Also we define the following quantities 

t^ln/' Rf" 

_ d\nf _ Rf _X4 
^~ d\nR~ f ~ X3' ^ ' 

Wc note that a constant value of m leads to the models with f{R) = ^1 + ^2R^~^"^ where the 
parameter m shows the deviation of the background dynamics from the standard model and ^1 
and ^2 are constants. However, in general the parameter m depends on R and R itself can be 
expressed in terms of the ratio r — ^ . This means that m is a function of r, that is, m = m(r). 
Based on the new variables, the Friedmann equation becomes a constraint equation so that we 
can express one of these variables in terms of the others. Introducing a new time variable 
r = In a = and eliminating xi (by using the Friedmann constraint equation) we obtain the 
following autonomous system 

— — = a;2(x5 + X4 + 2) , (16) 



dN 

dxs X4X5 
dN m 



+ x^{2xi + + 4) , (17) 



^ = ^ + ^4(2X4 + 4), (18) 

— = {X2 + X5){X5 + X4) + 1 - 3X3 - - 2X2 , (19) 

and 

xi = = 1 - X2 - X3 - X4 - X5 . (20) 

The deceleration parameter which is defined as q = —1 — ^ , now can be expressed as 

q=l+X4, (21) 
and the effective equation of state parameter of the system is defined by 

^eff = -1 - 1^ • (22) 
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Figure 1 : The effective equation of state parameter (left panel) and the deceleration parameter (right 
panel) of the critical point £ versus m. The shaded region denotes a non-phantom acceleration era 
which occurs in the parametric space with 0.45 < m < 0.67 . A slightly phantom-like behavior exists 
for m < 0.45 and when m — t- , it reaches to a strong phantom-like phase. For negative values of m 
the EoS parameter is stiff-like. 

2.3 Critical points and their stability 

The critical points of the scenario and some of their properties are hsted in table 1. In this 
table, F is defined as 

^_ 1 4m2-9m + 2± V-160 + 272 - 111 m2 + 4 m + 4 
^ 2 2 m2 - 3 m + 1 ' 

We consider only the plus sign of this equation in our forthcoming arguments. The minus sign 
does not create suitable cosmological behavior since it leads to Weff < — 10 or Weff > 0.7 for 
point £. 

In table 1, the critical points A , B and C are independent of the form of f{R)- Nevertheless, 
the stability of these points depends on the form of f{R) explicitly. The critical curve T> exists 
just for f{R) models with m(r = — |) = | (for instance, in models of the form f{R) = R+'jR~"' 
that m is defined as m(r) = _ "(^+'') ^ ^i^e critical curve V exists just for n = |). The value 
of the effective equation of state parameter corresponding to this critical curve depends on the 
coordinate X2 ■ Then, different intervals on this curve describe different era of the universe 
evolution (for example, a formal de Sitter-type era occurs for curve V at X2 = ~\)- Since the 
phase space behavior of the point S (with m(r) = 1 + r) depends on the form of f{R) via m, 
we treat it separately in section 3. In which follows, we classify the important subclasses in 
order to see their dynamical behaviors. 
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a) The radiation dominated era 

The points A and B demonstrate effectively the radiation dominated epoch of the universe 
{cueff = I with a{t) (X t^). This phase can be reahzed also by the curve V for = j ■ The 
stabihty of these points is investigated in which follows. 

b) The matter dominated era 

The matter era {uj^ff = 0) could be existed for models with m(r = — |) = | by the local- 
ized point = I on the curve V . Note that this matter era is properly described by a cosmic 

2 

expansion with scale factor a{t) = ao{t — to)^- This era also can be realized in localized point 
£ with m(r = -0.13) = 0.87. 

c) The de Sitter era 

The de Sitter phase (cUg// = —1) in the normal branch of this DGP-inspired f{R) model is 
realized by the curve C of critical points. It is important to point out that the mentioned de 
Sitter solution is the standard de Sitter phase just for X2 = 2 , since in this localized point the 
matter density parameter vanishes {Qm = 0). Also this phase can be realized from the curve 
of the non-localized points V with X2 = (this de Sitter point can be regarded as the special 
case of the curve C in the locahzed point X2 — — |). But one can see that this point gives no 
standard de Sitter era since its is non- vanishing (in this case — occurs at ^2 = |). 

d) Transition from g > to g < 

The critical point S depends explicitly on the form of f{R) via m and this is the case also 
for stability of this point. The point S describes a phase transition of the universe from decel- 
eration to the acceleration era at We// = — | and g = in this model for m(r = —0.33) = 0.67. 
Also the mentioned feature for f{R) models with m(r = ^|) = | (corresponding to the 
curve V) occurs at the fixed point x*2 = \ . So, one of the important feature of this model 
is that it clearly realizes the late-time acceleration of the universe in its normal branch for 
< m < 0.67 . In figure 1 the locahzed point S represents the deceleration phase for m > 0.67 
and a non-phantom accelerating phase for 0.45 < m < 0.67. For m < 0.45 , since g < — 1 , the 
model realizes an effective phantom phase with possibility of future big rip singularity which is 
characteristics of a non-canonical (phantom) field dominated universe. Similarly, the effective 
phantom behavior for the curve V occurs at 0^2 < — | . 

In the next step we determine the stability of the critical points under small perturbations. 
The stability of these points is determined by the eigenvalues of the Jacobian matrix. For a 
general f{R) term on the brane, stability of the critical points depends on the form of f{R) (or 
equivalently on the parameter m). It is obvious that in general m — m{r) is not a constant; 
it is a function of other variables so that one can expand this function of curvature about any 
of the fixed points. The results of our investigation for stability of critical points mentioned in 
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Table 1: Location, effective EoS and deceleration parameter of the critical points for the normal 
branch of a general DGP-inspired f{R) scenario. The fixed point T) exists only for those f{R) 
models that m(r — — ^- 
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2m 


r 
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table 1, are summarized as follows: 

• Point A 

As has been mentioned, this point is a radiation dominated era. In this case the eigenvalues 

are 

4(m-l) , , 

-2, ^, -4± i. (23) 

m 

Note that around this point, m = m_4 = m(r = 0) . Hence this point is a spiral attractor 
if < m{r = 0) < 1 , otherwise it is a saddle point. The corresponding 2D phase space for 
arbitrary m — m{r) is shown in figure 2 (left panel). Note that the Jacobian matrix in this 
point has no dependence on the m'(r) since this point lies around r — 0. Here a prime denotes 
derivative with the respect to r . 

• Point B 

This point is also a radiation dominated era in which the eigenvalues are as follows 



Ai,2 = --(11±\/T99z) 



3 — m + m'r(l + r) y ("^^ ~ "^)^ ^ m'rm'^{l — r) — 2 mm'r(l + r) + m' (r + r^)^ 

^3,4 = ?i ± 



m? rn? 



(24) 

where the parameter m should be expanded around the point B (that is, m = rrits). Here 
the stability issue depends also on m and m', so that this point can be either a saddle or a 
spiral attractor. The corresponding 2D phase space for arbitrary m = m{r) is shown in figure 
2 (right panel). An important point should be emphasized here: by setting X4 = in plotting 
figures 2, their dependence on m(r) is wasted. However, the situation is different if one plots 
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Figure 2: The phase subspace X3 - X5 of our setup at X2 = X4 = (left panel) and X2 = | and X4 = 
(right panel) . The critical point shown in the left panel is point A and in the right panel it is point 
B. Point A in the mentioned subspace is a spiral attractor and point ;B is a saddle point. Note that in 
the 4D phase space these are either spiral attractor or saddle point depending on the form of f{R) ■ 



the 3D subspace for these points. In this case, one should determine the form of the 

function m(r = ^). For a constant m, the eigenvalues reduce to (2 , — ± VT99 i) , '^^'^^'^^ ) 
which indicates that the point i3 is a saddle point for constant values of m. 



• Curve C 



Generally, if a nonlinear system has a critical curve, the Jacobian matrix of the linearized 
system at a critical point on the line has a zero eigenvalue with an associated eigenvector tangent 
to the critical curve at the chosen point. When dynamical variables are not independent, some 
eigenvalues of the Jacobian matrix are zero. In this case, the phase space of the nonlinear 
system reduces to a lower dimensional phase space. The stability of an specific critical point on 
the curve can be determined by the nonzero eigenvalues, because near this critical point there 
is essentially no dynamics along the critical curve (i.e., along the direction of the eigenvector 
associated with the zero eigenvalue). So, the dynamics near this critical point may be viewed 
in a reduced phase space obtained by suppressing the zero eigenvalue direction. On the other 
hand, such curves are actually normally hyperbolic [11,12]. We consider a point on the curve 
C with coordinates (2, 1, —2, 0). This point is a standard de Sitter phase. The eigenvalues 
corresponding to this point follows 



0. i(.-i^-4 

3^ X 



1/ 17 „\ , .v^/ 17x 
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X 



(25) 



where x is defined as 



X 



-460 m + 54 + 3 ^3^8019 - 1840 m + 108 ) 

m 



(26) 
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Figure 3: 2D phase space 3:3 - X5 for X2 = 1 and X4 = — 2 (left panel) and the 3D phase space - 
X3 - X5 for X4 = — 2 (right panel). The critical point which is shown in the left panel is a center and is 
related to the standard de Sitter phase. The critical points shown in the right panel which lie on the 
line X3 = |(11 — 4x2) are spiral attractors. These figures are plotted for m(r = —2) = \. 



The parameter ra should be expanded around the standard point of C as defined previously 
with mc = 'm{r = —2). Here the stability issue depends only on m. It is a stable spiral in the 
subspace of the last two eigenvalues when < m. On the other hand, the second eigenvalue is 
negative for m < and 0.09 < m. So the standard de Sitter point in the 4D phase space is a 
spiral attractor if 

m > 0.09 . (27) 

This point for other values of m is a saddle point. Figure 3 shows the 2D and 3D phase spaces 
of the critical curve C. We note that in the 2D subspace, the de Sitter curve is reduced to a de 
Sitter point (as figure 3 shows, this point is a center). Therefore, the center manifold theory 
is required to investigate its stability [13,14]. In 3D subspace (right panel), the non-localized 
points in the X2 - X3 plane lie on the line X3 = |(11 — 4x2). Figure 3 is plotted for m = | which 
satisfies (27), therefore this curve is a spiral de Sitter attractor. 



3 Analytical results for some specific models 

As we have mentioned previously, m is a function of r, that is, m = m{r). Since the ACDM 
model defined with f{R) = R — 2A, corresponds to m = 0, we can say that the quantity 
m characterizes the deviation of the background dynamics from the standard ACDM model. 
Now we consider two specific model of f{R) in order to obtain more obvious results. We also 
focus on the cosmological viability of these models. A cosmologically viable scenario contains 
an early time radiation dominated era followed by a matter dominated era that reaches to a 
standard de Sitter phase which is a stable attractor. We focus here only on the last two stages: 
matter domination and then a stable de Sitter attractor. At the first stage, the existence of 
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a matter dominated era {weff — 0) constrains our DGP-inspired f{R) model only to those 
f{R) functions that m(r = — |) = | and m{r = —0.13) = 0.87 (see table 1). In the first case, 
there exist two dc Sitter phases: the localized point X2 = — | on the curve V which is a formal 
de Sitter point, and also the de Sitter curve C which realizes a standard de Sitter point at the 
localized point X2 — 2. The second case is associated just to the de Sitter curve C. A correct 
connection between the unstable matter dominated era and the stable, standard de Sitter era is 
necessary condition for cosmological viability of a scenario. This connection can be investigated 
in the m — r plane. 



(30) 



A) f{R) = R + 7i?-" 

For this model, the parameter m takes the following form 

„„ = Z?1(1±I), (28) 
r 

which is independent on 7 . On the other hand, the point £ of table 1 is characterized by the 
following relation 

m(r) = r + 1 . (29) 

Equations (28) and (29) give two solutions nii — and 7712 = 1 — n for m. So, for the mentioned 
f{R) function, point £ is characterized by the following relation 

_ / r(^-) + 4(i-n) r(^-») + 4(i-n) r(^-) 

2n(l - n) ' 2(1^75 ' J' ^ ~^ ~ W^) 

where 1(1^") = r|„=i„„. 

Note that the critical point £ for m = (that is, £^(o)) is indefinite and therefore we exclude 
it from our considerations. Now we investigate the stability of the critical point £(i-n)- Since 
the EoS parameter corresponding to this point varies with m, in order to determine the stability 
of this point, one has to fix the value of m. 

In figure 4 (left panel), we have shown the behavior of the parameter m as a function of r for a 
special f{R) model given as /(i?) = R + 7/?" with n = 0.13 . As this figure shows, this model 
contains a connection between the matter era and the standard de Sitter era which is located 
at m = —0.065. However, this model is not cosmologically viable since it reaches an unstable 
standard de Sitter era. In figure 4 (right panel), we plotted the curve m(r) for f{R) — R + ^R"' 
with n = ^. The matter dominated era evolves to the standard de Sitter era at m = —0.25 
which indicates that the standard de Sitter era is unstable (see Eq. (27)). So, this model is 
not cosmologically viable too. Note that in the right panel of figure 4, the point A lies also on 
the dashed line which is corresponding to the critical point £. This feature indicates that this 
model contains the critical point £ with m(r — — ^ which is a non- phantom acceleration 
phase (compare this case with figure 1 that m — ^ lies in the shaded region). 
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Figure 4: The solid curves illustrate the diagram m — r of f{R) = R + 7/?"" with n = 0.13 (left 
panel) and with n = ^ (right panel). The dashed- line is corresponding to the point £ of table 1. The 
intersection point A shows a matter domination phase. The standard de Sitter phase is corresponding 
to the line r = —2 and dotted (blue) line indicates the line of stability. 



B) /(i?) = i?"exp(|) 

In this model the parameter m is defined as 

f \ n + r{2 + r) 

mir) = , (31) 

r 

which is independent on 77 . As has been pointed out previously, the matter dominated era can 
be achieved from critical point S in which the relative parameter is given by equation (29). 
Equating (29) and (31), we obtain two solutions for m as m = i±v^EM _ Now, the critical 
point S gives a matter dominated phase (with m{r = —0.13) = 0.87) if we set n = 0.35, 
that is to say, if f{R) = R exp(^) . In the model with f{R) = Rexp{^) (models with 
m{r = — i) = i), the matter dominated era is realized by the curve V and the critical point 
S plays the role of a non-phantom acceleration era. In figure 5 (left panel), the curve m(r) is 
plotted for f{R) = R exp(^) . This model reaches the standard de Sitter phase at m = 0.2 . 
Therefore, this point is stable since it belongs to the region defined by relation (27). Finally, 
we plot the m(r) curve for f{R) = i?exp(-^) as shown in the right panel of figure 5. In this case 
there is an acceptable connection between the matter dominated era and the standard de Sitter 
phase since the standard de Sitter point at m = | is a stable attractor. We note also that in 
this DGP-inspired f{R) model there is a non-phantom acceleration phase which emerges from 
the point S with m{r = — i) = |. This is corresponding to the point A of figure 5 (right panel). 



4 Summary and Conclusion 

In this paper we investigated cosmological dynamics of the normal DGP setup in a phase space 
approach where the induced gravity is modified in the spirit of /(i?)-theories. The motivation 
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Figure 5: The diagram m — r for f{R) = -R"exp(-^) with n = 0.35 (left panel). The dashed-line 
shows the m(r) function of the critical point £. The intersection point A shows a matter domination 
phase. This model reaches to the standard de Sitter phase at B with m = 0.2 which is a stable point. 
In this model the matter dominated era can be realized only by £. 

for this study within a dynamical system approach lies in the fact that recently it has been 
revealed that the normal, ghost-free DGP branch has the potential to explain late-time speed- 
up if we incorporate possible modification of the induced gravity in the spirit of /(-R)-theories. 
In this respect, a phase space analysis of the scenario would be interesting to reveal some 
aspects of this late-time behavior. Especially the stability of this late-time de Sitter phase is 
important to have a cosmologically viable solution. We applied the dynamical system analysis 
to achieve the stable solutions of the scenario in the normal DGP branch. We have shown 
that generally there are some fixed points that one of those is the standard de Sitter phase. 
Therefore, the normal branch of this DGP-inspired braneworld scenario realizes the late-time 
acceleration phase of the universe expansion. However, the stability of this point depends on the 
form of f{R) via the parameter m = -j^^- Then, we investigated the cosmological viability of 
these setups. A cosmologically viable scenario contains an early time radiation dominated era 
followed by a matter dominated era that reaches a standard de Sitter phase which is a stable 
attractor. Here we focused only on the last two stages: matter domination era followed by a 
stable de Sitter attractor. To be more specific, we considered two models with f{R) = R + jK^ 
and f{R) = -R"exp(-^) in our DGP-inspired setup. The condition for existence of the matter 
domination era restricted us to consider two cases n = 0.13 and n = | for the first model 
and two cases n = 0.35 and n = 1 for the second one. On the other hand, it is shown that 
the standard de Sitter phase is stable just for m > 0.09. So, since the first model reaches the 
standard de Sitter phase (which is determined by the line r = —2) at m < 0.09 , it is not a 
cosmologically viable model. However, since the second model reaches this phase at m > 0.09, 
it is a cosmologically viable model for n = 0.35 and n = 1. 
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